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On the integrability of the SU (N) Hubbard model 
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Abstract 

We exhibit explicitly the intertwiner operator for the monodromy matrices of the 
recent proposed SU(N) Hubbard model ||. This produces a new family of non-additive 
i?-matrices and generalizes an earlier result by Shastry H. 



M.J.Martins 

Instituut voor Theoretische Fysica, Universiteit van Amsterdam 
Valcknierstraat 65, 1018 XE Amsterdam, The Netherlands 

and 

Universidade Federal de Sao Carlos 

Departamento de Fisica 
CP. 676, 13560 Sao Carlos, Brazil 



For many years the one- dimensional Hubbard model has generated a great deal of attention, 
both due to the integrability properties Jl], |2|] and the promising role as a paradigm for strongly 
correlated electrons systems p, |4j . The one- dimensional Hubbard model has always appeared as 
an unusual and solitary example, outside from an integrable family of exactly solvable models. 
Very recently, however, Maassarani [Q proposed an integrable SU(N) generalization of the 
Hubbard model by using as building blocks the Boltzmann weights of a "free-fermion" SU(N) 
spin chain || . This author also constructed the Lax operator of the two-dimensional "covering" 
vertex model whose transfer matrix commutes with the SU(N) Hubbard Hamiltonian. The 
purpose of this note is to complement Maassarani's work, by presenting the quantum .R-matrix 
which intertwines two monodromy matrices of the SU(N) "covering" Hubbard model. This 
not only guarantees that transfer matrices with different spectral parameters commute but 
also paves the way for an algebraic Bethe Ansatz solution (see e.g. 0). In other words, we 
are interested in finding the f?-matrix which solves the Yang-Baxter algebra 

-Ri 2 (A, /i)£oi(A)£ 2(» = £ 2(»£oi(A).Ri2(A, p) (1) 

where Coi(X) is the SU(N) Lax operator proposed in ref. ||. 

For this purpose, we follow closely the strategy devised by Shastry || who solved similar 
problem in the special case of N = 2. First we apply a "gauge" transformation on the 
original SU(N) Lax operator proposed by Maassarani [Q, rewriting it as 

C»{\)=£&\\)£®(\)^&&° 1 * (2) 

In this equation £q^(A), k — 1,2 are two commuting copies of the vertex operator of 
the "free-fermion" SU(N) model |J. The interaction is via the "azimuthal" operator || 



= y^(e^ — e$ N ) while constraint h(X) satisfies the relation 
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where the constant U plays the role of an on-site Coulomb interaction 
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Here we would like to solve relation (1) and (2) for N > 3. As an example let us start 
with the case N = 3. For this model it is convenient to write the SU(3) "free-fermion" vertex 
operator as 
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where a(A) = cos(A) and 6(A) = sin(A). For this representation the interaction term reads 

/l 0\ 

=0-10 (5) 
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We recall that we have performed an appropriate canonical transformation on the original 
expressions for £ { £\X) and of ref. §. In this way the "free-fermion" operator becomes 
more symmetrical, helping us to better sort out the many different entries of equation (1). The 
next step is to make an Ansatz for -Ri2(A, p), guided by the results found by Shastry for N = 2 
[0. This suggests us to seek an expression for the f?-matrix having the following form 



R 12 (X, p) = 4J(A - /x)££'(A - /i) + /(A, p)£%(X + p)p\ L, C^{X + p)p\ 
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The coupling f(X,p) is determined via a brute force analysis, i.e. by inserting the above 
Ansatz and the expression for the Lax operator in the the Yang-Baxter algebra (1). Here one 
can use useful identities between the "free-fermion" Lax operator, and for instance we have 



4?(A - p)£\l\X)prc%(p)p?> = C^>(p)p^L\l>(X)prc\l>(X - p) 
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Now we can put all these information together, by expanding the exponential part of the 
Lax operator (2), and by considering the above identities and its variants. After some algebraic 
manipulations we find that a necessary condition for this function is 

/(A, n) sin(A + fi) [tanh[/i(A)] tanh[/*(^)] + 1] = sin(A - n) [tanh[/i(A)] + tanh[/i(^)]] (10) 

Keeping in mind the constraint (3), it is not difficult to find that the above equation is 
satisfied by the following expression 



Remarkably enough, the solution for function /(A, fi) is precisely the one found by Shastry 
for N = 2. This means that the functional expression for the i?-matrix is indeed universal, and 
the structure of the Boltzmann weights depends only on a given representation for the "free- 
fermion" building blocks. To verify this fact we have performed extensive checks for N = 3,4 
and we found that solution (11) is also a sufficient condition. This leads us to conjecture that 
this should be valid for arbitrary integer N > 2. All this resembles much the Baxterization 
procedure, and it seems that a better understanding of these results should necessarily start 
with the discovery of the algebraic structure, such as Temperley-Lieb and braid-monoid, which 
is behind the i?-matrix construction (6) and (11). This could also prompt new representations 
for the building blocks, leading us to generate many more integrable "Hubbard models " . 
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In summary we have exhibit the f?-matrix for the recent proposed SU(N) Hubbard model. 
We also verified, for N = 3,4, that the f?-matrix solution satisfies the Yang-Baxter equa- 
tion, by following the approach of Shiroishi and Wadati ||. We find that the amplitudes 
for the "tetrahedral" Zamolodchikov algebra is precisely the same of that found for N = 2 
0, corroborating the universal scenario proposed above. We hope that such i?-matrix solution 
can be used together with the framework developed in ref. M to solve the SU(N) Hubbard 
model via an algebraic Bethe Ansatz method. 
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